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Abstract 

We study the correspondence between the interacting viscous ghost dark energy model 
with the tachyon, K-essence and dilaton scalar field models in the framework of Einstein 
gravity. We consider a spatially non-flat FRW universe filled with interacting viscous ghost 
dark energy and dark matter. We reconstruct both the dynamics and potential of these scalar 
field models according to the evolutionary behavior of the interacting viscous ghost dark 
energy model, which can describe the accelerated expansion of the universe. Our numerical 
results show that the interaction and viscosity have opposite effects on the evolutionary 
properties of the ghost scalar filed models. 
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1 GDE scenario 

The GDE density is proportional to the Hubble parameter [1, 2, 3, 4, 5, 6, 7, 8, 9] 

p D = aH, (1) 

where a is a constant with dimension [energy] 3 , and roughly of order of Aq CD where Aqcd ~ 
lOOMeV. 

Here we consider a spatially non-flat FRW universe filled with GDE and DM. Within the 
framework of FRW cosmology, the first Friedmann equation takes the form 

H2 + ^ = if,(Pn+Pm), (2) 

where M p = (87rG) -1 / 2 is the reduced Planck mass. Here k = 0, 1, —1 represent a flat, closed 
and open FRW universe, respectively. Also pD and p m are the energy densities of GDE and 
DM, respectively. 

Using the dimensionless energy densities defined as 

O _ Pm _ Pm _ PD_ _ PD n _ k 

m ~ Pcr " 3M2JJ2 ' p CI 3M*H* ' ^ " a 2 ff 2 ' ^ J 

the Friedmann equation (2) can be rewritten as 

i + n k = n D + n m . (4) 

Substituting Eq. (1) into p D = 3M£H 2 tt D yields 



n ° = sMh- (5) 



a 

Using the above relation, the curvature energy density parameter can be obtained as 



^/«fi*W = ^W, (6) 



v a 2 I \ a J \^£>o/ v u / 

where we take ao = 1 for the present value of the scale factor. 

Here, we extend our study to the viscous model of GDE. In the presence of viscosity, the 
effective pressure of DE takes the form 

p D = PD - (7) 

where £ = eH~ 1 pD is the bulk viscosity coefficient in which e is a constant parameter [10]. A 
viscosity e > will be able to drive acceleration [10]. 

We further assume the viscous GDE interact with DM. The interaction between DE and DM 
can be detected in the formation of LSS like the galaxy cluster Abell A586 [11]. In the presence 
of interaction, the continuity equations are 

Pd + 3H{1+lu d )pd = 9eH PD -Q, (8) 
p m + 3Hp m = Q, (9) 
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where up = Pd/pd is the equation of state (EoS) parameter of the interacting viscous GDE 
and Q stands for the interaction term. Following [12], we shall assume Q = 3b 2 H(p m + po) with 
the coupling constant b 2 . 

Taking time derivative of Eq. (1) and using Eqs. (2), (4), (5) and (9) gives 



PD 



n D -i-^ + b 2 (i + n k ) 

2-Q D 



(10) 



Taking time derivative of Eq. (5) and using (1) and (10) one can obtain the evolution of the 
interacting viscous GDE density parameter as 



d In a \Qd 



n D -i-^ + b 2 (i + n k ) 



(ii) 



which is same as that obtained for the interacting GDE in non-flat universe in the absence of 
viscosity [6]. It is interesting to note that the viscosity constant e does not affect the evolution 
of the GDE density parameter (11). Substituting Eq. (6) into (11) yields a differential equation 
for Qd(cl) which can be solved numerically with a suitable initial condition like Qd = 0.72. The 
numerical results obtained for Qr>(a) are displayed in Fig. 1 for different coupling constant b 2 . 
Figure shows that: i) for a given b 2 , increases when the scale factor increases, ii) At early 
and late times, increases and decreases with increasing b 2 , respectively. 

Substituting Eq. (10) into (8) gives the EoS parameter of the interacting viscous GDE model 

as 

1_ % + 262/1+^ 



= n D -2 +3e ' (12) 

which shows that in the absence of interaction and viscous terms, i.e. b 2 = e = 0, at early 
(Qd —> 0) and late (Qd 1) times ojd goes —1/2 and —1, respectively, and cannot cross the 
phantom divide line [4]. For the present time (do = 1), taking Qd = 0.72 and $\ = 0.01 [13] 
Eq. (12) gives 

co Dq = _o.78-2.196 2 + 3e, (13) 

which clears that for e = the phantom EoS parameter (ojd < — 1) can be obtained provided 
b 2 > 0.1. This value for b 2 is consistent with the recent observations in which we have that 
b 2 could be as large as 0.2 [14]. Also the phantom divide crossing is in accordance with the 
observations [15]. 

The evolution of the EoS parameter (12) for different b 2 and e is plotted in Figs. 2 and 3, 
respectively. Figure 2 shows that: i) for b 2 = 0, ojd decreases from —0.5 at early times and 
approaches to —1 at late times, ii) For b 2 / 0, ujd increases at early times and decreases at 
late times. The results of ujd in the absence of viscosity (e = 0) are in agreement with those 
obtained by [6]. Figure 3 clears that: i) for a given e, ujjj decreases with increasing the scale 
factor, ii) For a given scale factor, ujd increases when e increases. 

2 Ghost tachyon model 

The tachyon field is another approach for explaining DE. The tachyon energy density and pres- 
sure are [16] 

PT = ^L, (14) 
1-6 2 



3 



The tachyon EoS parameter yields 



PT = -V(<P)Jl-^. 



lot = — = 9 ~ 1- 



(15) 



(16) 



To reconstruct the tachyon filed via the interacting viscous GDE, equating (12) with (16), 
i.e. wd = wy, gives 

l_ik + 26 2(l+^ 



+ 3e = 2 - 1 



fi D -2 

Also equating Eq. (1) with (14), i.e. po = pt, gets 

V(<f>) 



aH 



1-6 2 



(17) 



(18) 



From Eqs. (17) and (18), the kinetic energy and potential of the tachyon field can be obtained 
as follows 

^-l-% + 26 2 (±^) 



J.2 



v{4>) 



or 



n D - 2 



1-^ + 26 



+ 3e, 



2( i±n 
n 



2-n D 



3e 



1/2 



(19) 



(20) 



Note that Eqs. (19) and (20) for the flat case, i.e. = 0, and in the absence of viscosity (e = 0) 
reduce to the results obtained by [8]. 

From Eq. (19) and using (5), one can get the evolutionary form of the ghost tachyon scalar 
field as 



vi) = — / n D 

a Ji 



ft D -l-% + 2o 2 



2/ i+n k 



Q D -2 



+ 3e 



1/2 



da 
a 



(21) 



where we take ao = 1 for the present time. The evolution of the ghost tachyon scalar filed, Eq. 
(21), for different values of b 2 and e is plotted in Figs. 4 and 5, respectively. Figures clear that: 
i) for a given b 2 or e, (j)(a) increases with increasing the scale factor, ii) For a given scale factor, 
4>{a) decreases and increases with increasing b 2 and e, respectively. Note that Fig. 4 shows only 
the real scalar field, i.e. 4> 2 > 0. Indeed, for b 2 = 0, 0.02 and 0.04 the scalar field 4> becomes 
pure imaginary ((f> 2 < 0) at a > 43.5, 2.8 and 2.1, respectively, and it does not show itself in Fig. 
4. To investigate this problem in ample detail, the evolution of the ghost tachyon kinetic energy 
X = 4> 2 /2, Eq. (19), for different values of b 2 and e is plotted in Figs. 6 and 7, respectively. 
Figure 6 confirms that for b 2 = 0, 0.02 and 0.04 the kinetic energy becomes negative (x < 0) at 
a > 43.5, 2.8 and 2.1, respectively. Figures 6 and 7 show that: i) for a given b 2 or e, the kinetic 
energy x decreases when the scale factor increases, ii) For a given scale factor, the kinetic energy 
decreases and increases with increasing b 2 and e, respectively. 

It is worth to note that from Eq. (14) due to having a real tachyon energy density we need to 
have (p 2 < 1 which is in accordance with Figs. 6 and 7. Moreover, from Eq. (16) for cp 2 < and 
< 4> 2 < 1 we have lot < — 1 and —1 < cot < 0, respectively, corresponding to the phantom [17] 
and quintessence [18] DE, respectively. In the absence of interaction (b 2 = 0), the kinetic energy 
of the ghost tachyon scalar field is always positive (see Fig. 7) and behaves like quintessence DE 
with lot = > — 1 (see Fig. 3). 
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The ghost tachyon potential, Eq. (20), versus the scalar field (21) for different b 2 and e 
is plotted in Figs. 8 and 9, respectively. Figures illustrate that: i) for a given b 2 or e, V{(j>) 
decreases with increasing <j>. This behavior is in agreement with the scaling solution V(4>) oc <p~ 2 
obtained for the tachyon filed corresponding to the power law expansion [19]. ii) For a given 
scalar field, V{4>) increases and decreases with increasing b 2 and e, respectively. 



3 Ghost K-essence model 

The K-essence scalar field model of DE is characterized by a scalar field with a non-canonical 
kinetic energy. The most general scalar field action which is a function of <f> and x = 4> 2 is 
given by [20, 21] 

5 = | d 4 x^p(^ X ), (22) 

where the Lagrangian density p(4>, x) corresponds to a pressure with non-canonical kinetic terms 
as 

P(^x) = fm-X + X 2 ), (23) 
and the energy density of the K-essence field <fi is 

P^,x) = fm-X + 3x 2 )- (24) 
The K-essence EoS parameter takes the form 

PW,X) 3x-l 
Equating (25) with (12), ujk = ujd, we get 

3 - % + 2b 2 (^pA -n D + 3e(n D - 2) 

x = v > . (26) 

5 - n k + 66 2 (i±^) - n D + 9e(fl D - 2) 



Using Eq. (26) and <p 2 = 2 X , we obtain the ghost K-essence scalar field as 

4>{a) - 0(1) = — / flf) 

a Ji 



5-n k + 6b 2 (^) -fl D + 9e(n D -2) 



which its evolution for different b 2 and e is displayed in Figs. 10 and 11, respectively. Figures 
present that: i) for a given b 2 or e, (f>(a) increases with increasing the scale factor, ii) For a given 
scale factor, <p(a) decreases and increases with increasing b 2 and e, respectively. 

The evolution of the ghost K-essence kinetic energy, Eq. (26), for different values of b 2 and 
e is plotted in Figs. 12 and 13, respectively. Figures clarify that: i) for a given b 2 or e, the 
ghost K-essence kinetic energy like the tachyon filed decreases when the scale factor increases, 
ii) For a given scale factor, the kinetic energy of the ghost K-essence filed like the tachyon model 
decreases and increases with increasing b 2 and e, respectively. If we compare Fig. 12 with 6 
we see that the kinetic energy of the ghost K-csscncc model in contrast with the ghost tachyon 
field is always positive. Note that the result of Fig. 12 is in contrast with that obtained by [9] 
who showed that for a given b 2 , the kinetic energy of the ghost K-essence filed increases with 
increasing the scale factor. This difference may come back to this fact that the K-essence filed 
selected by [9] is a purely kinetic model in which the action (22) is independent of 4>. This yields 
the energy density and pressure of a purely kinetic K-essence which are different from those 
considered in Eqs. (23) and (24). 
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4 Ghost dilaton model 



The dilaton scalar field model is also an interesting attempt to explain the origin of DE using 
string theory. The dilaton pressure and energy density are given by 



p D = -x + ce A< V, 



PD 



= - X + 3ce A V 



(28) 
(29) 



where c and A are positive constants and % 



WD 



Pd_ 

PD 



b 2 /2. The dilaton EoS parameter takes the form 
1 



(30) 



Equating (30) with (12) gives the solution 

3 



ce X(t> x 



% + 26 : 



2 1 i+n k 
n D 



fl D + 3e(n D - 2) 



5 - O fc + 66 2 (i±g^) -n D + 9e(n D -2) 



(31) 



then with the help of x = <f/2, we obtain 



e 2 ( 



3 - % + 26 2 (i±gt) - n D + 3e(n D - 2) 
5-ft k + 66 2 (±±gt) - n D + 9e(n D - 2) 



1/2 



(32) 



Finally we obtain 
4>{o) = -r In ( 



A 



e 2 + 



3M p 2 A 
2ay / c 



2Q t 



+ 46 5 



2^ D + 6e(fi D - 2) 



5 _ n k + 66 2 - tt D + 9e(n D - 2) 



1/2 



da 
a 



> . 



(33) 

The evolution of the ghost dilaton scalar field (33) for different b 2 and e is displayed in Figs. 
14 and 15, respectively. Figures present that: i) for a given b 2 or e, 4>(a) increases with increasing 
the scale factor, ii) For a given scale factor, 4>(a) decreases and increases with increasing b 2 and 
e, respectively. 

With the help of Eq. (31) we plot the evolution of the ghost dilaton kinetic energy for 
different b 2 and e in Figs. 16 and 17, respectively. Figures show that for a given b 2 or e, the 
kinetic energy of the ghost dilaton field like the tachyon and K-essence models decreases with 
increasing the scale factor. 



5 Conclusions 

Here we investigated the interacting viscous GDE model in the framework of standard FRW 
cosmology. We considered a spatially non-flat FRW universe filled with GDE and DM. We 
obtained the evolution of the fractional energy density and EoS parameters of the interacting 
viscous GDE model throughout history of the universe. Furthermore, we reconstructed both the 
dynamics and potential of the tachyon, K-essence and dilaton scalar filed models of DE according 
the evolutionary behavior of the interacting viscous GDE model. Our numerical results show 
the following. 
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(i) The evolution of the interacting viscous GDE density parameter Qjj is independent of 
viscosity constant e. But for a given coupling constant b 2 ,Q.D increases with increasing the scale 
factor. Also at early and late times, Sip increases and decreases, respectively, with increasing 
b 2 . 

(ii) The EoS parameter ojjj of the GDE model in the absence of viscosity, can cross the 
phantom divide line (00 d < —1) at the present provided b 2 > 0.1 which is compatible with the 
observations. Also in the absence of viscosity for a given coupling constant b 2 , u>z> increases and 
decreases at early and late times, respectively. Moreover, in the absence of interaction for a 
given viscosity constant e, ojd decreases when the scale factor increases. For a given scale factor, 
wd increases with increasing e. 

(iii) The ghost tachyon scalar filed for a given b 2 or e, increases with increasing the scale factor. 
Also for a given scale factor, it decreases and increases with increasing b 2 and e, respectively. 
For a given b 2 or e, the ghost tachyon kinetic energy %(a) and potential V{4>) decrease with 
increasing the scale factor and scalar filed, respectively. For a given scale factor, %(a) decreases 
and increases with increasing b 2 and e, respectively. For a given scalar field, V{4>) increases and 
decreases with increasing b 2 and e, respectively. 

(iv) The ghost K-essence scalar filed for a given b 2 or e increases with increasing the scale 
factor. But its kinetic energy decreases. For a given scale factor, the ghost K-essence scalar filed 
decreases and increases with increasing b 2 and e, respectively. This behavior also holds for the 
kinetic energy of the ghost K-essence model. 

(v) The ghost dilaton scalar filed and its corresponding kinetic energy for a given b 2 or e 
behave like the ghost K-essence model. 

All mentioned in above illustrate that the interaction and viscosity have opposite effects on 
the dynamics of ghost tachyon, K-essence and dilaton scalar field models of DE. 
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Interacting viscous GDE model 



1.0 




0.8 






0.6 






0.4 




=0.0 




=0.02 






=0.04 ; 


0.2 






0.0 







1 2 3 4 5 



a 

Figure 1: The evolution of the GDE density parameter, Eq. (11), for different coupling constants 
b 2 . Auxiliary parameters are Qd = 0.72 and ftk = 0.01. 



Interacting GDE model (e = 0) 




12 3 4 5 



Figure 2: The evolution of the EoS parameter of GDE, Eq. (12), for different coupling constants 
b 2 with e = 0. Auxiliary parameters as in Fig. 1. 



Viscous GDE model (b 2 = 0) 
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Figure 3: Same as Fig. 2 for different viscosity constants e with b 2 = 0. Auxiliary parameters 
as in Fig. 1. 
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Ghost tachyon model (e = 0) 
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Figure 4: The evolution of the ghost tachyon scalar filed, Eq. (21), for different coupling 
constants b 2 with e = 0. Auxiliary parameters are Qr> = 0.72, f^ = 0.01 and 0(1) = 0. 
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Figure 5: Same as Fig. 4 for different viscosity constants e with b 2 = 0. Auxiliary parameters 
as in Fig. 4. 
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Ghost tachyon model (e - 0) 
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Figure 6: The evolution of the ghost tachyon kinetic energy \ = 2 /2, Eq. (19), for different 
coupling constants b 2 with e = 0. Auxiliary parameters as in Fig. 4. 



Ghost tachyon model (b 2 = 0) 
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Figure 7: Same as Fig. 6 for different viscosity constants e with b 2 = 0. Auxiliary parameters 
as in Fig. 4. 



11 




Figure 8: The ghost tachyon potential, Eq. (20), versus the scalar field <j> for different coupling 
constants b 2 with e = 0. Auxiliary parameters as in Fig. 4. 



Ghost tachyon model (b 2 = 0) 




Figure 9: Same as Fig. 8 for different viscosity constants e with b 2 = 0. Auxiliary parameters 
as in Fig. 4. 
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Ghost K-essence model (e = Q) 
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Figure 10: The evolution of the ghost K-essence scalar field, Eq. (27), for different coupling 
constants b 2 with e = 0. Auxiliary parameters as in Fig. 4. 



Ghost K-essence model (b 2 = 0) 
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Figure 11: Same as Fig. 10 for different viscosity constants e with b 2 = 0. Auxiliary parameters 
as in Fig. 4. 
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Ghost K-essence model (e = 0) 
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Figure 12: The evolution of the ghost K-essence kinetic energy x = 2 /2, Eq. (26), for different 
coupling constants b 2 with e = 0. Auxiliary parameters as in Fig. 4. 



Ghost K-essence model (b 2 = 0) 
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Figure 13: Same as Fig. 12 for different viscosity constants e with b 2 = 0. Auxiliary parameters 
as in Fig. 4. 
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Ghost dilaton model (e = 0) 
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Figure 14: The evolution of the ghost dilaton scalar field, Eq. (33), for different coupling 

constants b 2 with e = 0. Auxiliary parameters are Qd = 0.72, Qk = 0.01, (f)(1) = and 
3M p x _ 

2a y/c ~ 



Ghost dilaton model (b 2 = 0) 
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Figure 15: Same as Fig. 14 for different viscosity constants e with b 2 = 0. Auxiliary parameters 
as in Fig. 14. 
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Ghost dilalon model (e = 0) 
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Figure 16: The evolution of the ghost dilaton kinetic energy \ = V ! /2, Eq. (31), for different 
coupling constants b 2 with e = 0. Auxiliary parameters as in Fig. 14. 



Ghost dilaton model (b 2 = 0) 



0.7 






0.6 


\ • 

- \ \ 




0.5 


: \\ \ 
" v • 






■ V *. 

■ \\ • 


— e=0.0 


0.4 


- V ** 
V \ 
■ v * 


- - e=0.1 . 
- e=0.2 : 


0.3 


\> ** 




0.2 






0.1 







1 2 3 4 5 



a 

Figure 17: Same as Fig. 16 for different viscosity constants e with b 2 = 0. Auxiliary parameters 
as in Fig. 14. 
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